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ABSTRACT 

We investigate the shapes and mutual alignment of voids in the large scale matter 
distribution of a ACDM cosmology simulation. The voids are identified using the novel 
WVF void finder technique. The identified voids are quite nonspherical and slightly 
prolate, with axis ratios in the order of c : b : a w 0.5 : 0.7 : f . Their orientations are 
strongly correlated with significant alignments spanning scales > 30/i -1 Mpc. 

We also find an intimate link between the cosmic tidal field and the void orien- 
tations. Over a very wide range of scales we find a coherent and strong alignment of 
the voids with the tidal field computed from the smoothed density distribution. This 
orientation-tide alignment remains significant on scales exceeding twice the typical 
void size, which shows that the long range external field is responsible for the align- 
ment of the voids. This confirms the view that the large scale tidal force field is the 
main agent for the large scale spatial organization of the Cosmic Web. 

Key words: Cosmology: theory - large-scale structure of Universe - Methods: data 
analysis - numerical 



1 INTRODUCTION 

Galaxy redshift surveys and computer simulations of cos- 
mic structure formation have sho wn that mat t er tra ces out 
a frothy pattern, the Cosmic Web. lBond et al\ <| 19961 ) estab- 
lished how the tidal forces induced by the inhomogeneous 
cosmic matter distribution are the main agent shaping the 
Cosmic Web. The resulting cosmic tidal force field is the 
source of the large scale, coherent, anisotropic forces which 
shape the cosmic matter distribution into characteristic fil- 
amentary and sheetlike patterns. 

Voids are a dominant componen t of 



the Cosmic Web (see e. g . iTullv et al. I 120071 : 
iRomano-Dfaz fc van de Weygaertl I2007D . occupying most 
of the volume of space. In this paper we present evidence 
for significant alignments between neigbouring voids and 
establish the intimate dynamic link between voids and 
the cosmic tidal f orce field. The Watershed Void Finder 
(WVF) procedure (|Platen. van de Weygaert fc Jonesll2007l ') 
that we use to identify voids is crucial to the succes of this 
analysis. The WVF technique is based on the topological 
characteristics of the spatial density field and thereby 
provides objectively defined measures for the size, shape 
and orientation of void patches. 

Large voids form around deep density troughs in the 
primordial density field. The main aspects of the evolu- 
tion of a large void may be understood on the basis of the 
expansion of simple spherical and isolated under-densities 
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(e.g. lBertschingerl (|l985h V Under-dense regions expand with 
respect to the background Universe and in general will 
have the tendenc y to grow more s pheric al wit h time (jlckd 
Il984h However, IShandarin et all (|2006t ) and iPark fc Led 
(2007a) have emphasized that, in realistic cosmological cir- 
cumstances, voids will be nonspherical. This is quite ap- 
pare nt in images of. for example, the Millennium simula- 
tion ijSpringel et alj|2005l ). Moreover, substructures within 
voids display a manifest alignment along a preferred direc- 
tion. This is, in part, a consequence of the relatively strong 
influence of the surrounding inhomogeneous matter distri- 
bution on the void's structure and evolution. 

Voids evolve hierarchically: as they expand with 
respect to the background Universe they merge with 
their surroundin g peers , building voids of ever larger sizes 
( Rcgos fc Celled Il99ll ; Ivan de Weygaert fc van KamperJ 



119951 : iGottlober et al.l 120031 ; IColberg et all 120051 ). Small 
voids can survive as substructure within the larger 
voids or may disappear under the gr avitational impact 
of surrounding overdense stru ctures (|Sahni et alj 1 19941 : 
ISheth fc van de Weygaert] |2004| ). The result is a void pop- 
ulation whose scale is distributed around a characteristic 
void size: we shall show this using our WVF sample. 

A major manifestation of large scale tidal influences is 
that of the alignment of shape and angular m omentum of 
objects (see lBond et <aZ.lll996l ; [Pes iacqucs 2 0071) . The align- 
ment of the orientations of galaxy haloes, galaxy spins and 
clusters with larger scale structures such as clusters, fila- 
ments and superclu sters have been the subject of numer- 
ous studies (see e.g. iBinggcli 1982]; lBondlll987l ; iRhee et all 
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Figure 1. The four panels show a slice through the density field with the void-segmentation superposed. The first panel (a) a slice through 
the entire simulation. The image shows the density field of the simulation. The density field is determined by DTFE. Superimposed are 
the black solid lines indicating the boundaries of the WVF voids. The white box indicates the section which forms the subject of the 
next three frames. Second panel: density field in the subbox with the WVF void boundaries superimposed (solid lines). The blue bars 
indicate the direction of the projection of the longest void axis. Third panel: the WVF void boundary contours (solid lines) with the fitted 
void shape ellipsoids sperimposed. Fourth panel: on the landscape with WVF void boundaries the tidal field compressional component 
is represented by tidal bars (red), representing the direction and strength of the tidal field. Also depicted are the void shape bars (blue). 



199ll: IPlionis fc Basilakod 120021: iBasilakos et all l200fj; 
Truiillo et al.ll2006l:lAragon-Calvo et al.ll2007l ; lLee fc Evrardl 
20071 : iPark fc Ledl2007bl ; ILee et al.ll2007M. Re cent analytical 
and numerical work by IPark fc Led (|2007ah discussed the 
ma gnitude of the tidal contribution to the shape of voids (see 
also ILee fc ParkboOot ) . Describing the evolution of voids un- 
der the influence of tidal fields by means of the Lagrangian 
Zel'dovich approximation, they found that the ellipticity dis- 
tribution of voids is a sensitive function of various cosmo- 
logical parameters and remarked that the shape evolution of 



voids provides a remarka bly robust constr aint on the dark 
energy equation of state (|Lee fc Parkll2007h . 

In section 2 we introduce the N-body simulation which 
we used to identify the voids, along with a brief descrip- 
tion of the WVF method we use to identify the voids in 
our sample. The results for individual void characteristics, 
in particularly the size and shape of voids, are presented in 
section 3. Then in section 4 we present the analysis of the 
alignment of voids in the sample. The possible origins for 
the void alignment are discussed in section 5, and in section 
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6 we focus on the role of tidal fields in determining the evo- 
lutionary history of voids. Finally, section 7 concludes with 
a short discussion and prospects. 



2 THE VOID SAMPLE 

It is essential to have a statistically representative void sam- 
ple for our study of the shapes of voids and their relative 
alignment. We also need a sufficiently large spatial cover- 
age in order to evaluate tidal field contributions over a large 
range of scales. Thus it is necessary to have a simulation 
with a large box size. We used the VLS simulation from the 
Virgo Supercomputing Consortium 1 , a simulation of cosmic 
structure formation in a ACDM cosmology with fl m ,o = 0.3, 
^a,o = 0.7, H = 70kms~ 1 Mpc- 1 and a a = 0.9. The sim- 
ulation contains 512 3 N-body particles, each with a mass 
of 6.86 x I0 10 h~ x Mq in a comparatively large box of side 
479/i _1 Mpc. We used the particle distribution at the final 
timestep, corresponding to the current cosmic epoch a = 1. 

2.1 Void Identification 

The particle distribution of the simulation is converted into a 
grid-based density fi eld by means of Delaunay Tr iangle Field 
Estimator (DTFE) (|van de Weygaert fc Schaadl2007l ). The 
key feature of DTFE is that it is a reconstruction method 
that does not involve any user-defined parameters. More- 
over, its natural and intrinsic filtering and interpolation pro- 
cedure yields a density field that remains faithful to both the 
density and shape structure of the original particle distri- 
bution. DTFE guarantees an optimal representation of the 
geometry and topology of the weblike patterns in the mass 
distribution. A slice through the DTFE resampled density 
field is shown in the first panel of figure Q] 

After resampling the point density distribution with 
DTFE, we applied the Watershed Void Finder (WVF) 

1 The Virgo Consortium VLS ACDM data is available at 
http: / / www.mpa-garching. mpg.de / galform / virgo /vis / 



l|Platen. van de Weygaert fc Jonesl I2007T ) in order to iden- 
tify the voids in the simulation. The WVF procedure traces 
the basin and its encompassing boundary for each minimum 
in the spatial mass distribution. Because WVF is free of any 
a priori shape criterion it is ideally suited to analyzing the 
morphology of underdense regions in the cosmic web. 

The WVF procedure w o rks as follows (see 
IPlaten. van de Weygaert fc Jonesl 120071 . for a detailed 
description). First we remove shot noise effects in the raw 
DTFE density field by smoothing at the Nyquist scale 
(~ l.O/i, Mpc); this is done using a Butterworth n — 2 
filter. Then we identify the minima in the resulting density 
field and tag each minimum with a unique identification 
number. We can visualize the density field as a landscape 
with basins around each minimum. Starting from the 
minima we flood this landscape. As the flood level rises the 
basins surrounding each miminum start to fill up. Walls 
gradually emerge along the crest lines of the density field at 
the interstices between distinct regions. The final product 
of the procedure will be a division of the density field into 
individual void cells. 



3 VOID CHARACTERISTICS 

In the first panel of figure [TJ each void identified by WVF 
is indicated by the solid black line marking the surround- 
ing void boundary. The image shows that we indeed have 
a large number of voids in our sample. In the upper right- 
hand panel we zoom in on a region of the simulation box (the 
white box in first panel) which we use for the purpose of il- 
lustrating the subsequent steps (the analysis is done on the 
entire box). The upper right panel also has lines indicating 
the (projected) magnitude and direction of the major axis of 
each void. Visual inspection of the void regions in this panel 
manifestly suggests some degree of alignment. Quantitative 
analysis will establish the significance of that alignment. 
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3.1 Void Size 

The size and shape distribution of the WVF-identified 
voids is shown in figure [2] In the simulation we found 
29000 voids, with an average void size Rv ~ 8.5/i Mpc. 
The void size distribution in figure [2^, clearly shows a 
peak in the size distribution, co nfirming the prediction by 
ISheth fc van de Weveaer tl l|2004l) . 

3.2 Void Shape 

For each individual void region we calculate the shape-tensor 
Sij by summing over the N volume elements k located 
within the void, 

Sij =- ^2x ki x kj (offdiagonal) (I) 

S " = ( x * _ x **) (diagonal) , 

fc 

where Xfc is the position of the fe-th volume element within 
the void, with respect to the (volume-weighted) void center 
r^j i.e. Xfc = rfc v v . The shape tensor Sij is related to the in- 
ertia tensor Ty. However, it differs in assigning equal weight 
to each volume element within the void region. Instead of bi- 
asing the measure towards the mass concentrations near the 
edge of voids, the shape tensor Sij yields a truer reflection 
of the void's interior shape. 

When evaluating alignments we will mostly use the 
traceless version of the shape tensor: 

Sij = Sij — —SSij , S = ^ ^ iSfcfc (2) 

fc 

The length and orientation of the shape ellipsoid are inferred 
from the eigenvalues and eigenvectors of Sij . The eigenvalues 
of the void ellipsoid are labelled so that si > S2 > S3. The 
corresponding semi-axes are oriented along the directions of 
the corresponding unit eigenvectors b s i, e S 2 and e„3. The 
length of the semi-axes a; are 

2 5 

a i = 2N { S2 + S3 ~ Sl J ( 3 ) 

with cyclic permutation for the other axes (N is the number 
of volume elements within the void). The longest axis is 
directed along e s i, the shortest along e s ^. 

We quantify the shape of the void ellipsoids in terms of 
the two axis ratios 7721 = 02/01 and 7732 = 03/02. Amongst 
the triaxial void ellpsoids we make a distinction between 
oblate ones, with 7)21 > 7)32, and prolate ones having 7721 < 
7732- If voids were spherical we would have 7721 = 7732 = 1 
(though, in reality, perfect sphericity does not exist). 

3.3 Results on Void Shapes 

In the density field in the upper righthand frame of figure [T] 
we have indicated the projected direction of the longest axis 
of the voids by means of a bar along the corresponding eigen- 
vector, with the length of the bar proportional to the size 
of the void axis. It allows a superficial comparison with the 
related WVF void patches (indicated by the solid black lines 
marking their boundary). 

The lower lefthand panel of figure [1] shows the shape- 
ellipsoids of each void within the region of the simulation 



box shown in the top-right panel. The centers of the shape- 
ellipsoids are located at the center of each void with their size 
scaled according to match the void volume: they give a con- 
venient visual impression of the void shapes and sizes. The 
reasonably accurate degree to which the ellipsoids reflect the 
shape of the voids may be inferred from a comparison with 
the WVF void regions themselves. 

Clearly, the voids are quite nonspherical. A more quan- 
titative impression of the shape distribution may be found 
in the second and third panel of fig. [2] which show that the 
voids are far from spherical. 

The intrinsic triaxial shape of voids is apparent from 
the distribution in the second panel of fig. [2] This shows a 
scatter plot of of the two axis ratios, 7712 versus 7723 for every 
void in the sample. To guide the eye we have superimposed 
the inferred isodensity contours defined by the point density 
in the scatter plot. We find a slight asymmetry: there are 
more prolate voids than oblate ones. 

The third panel shows that the distribution of the el- 
lipticity e = 1 — 0,3 /a-i is skewed towards higher values of 
e > 0.5. The void population is marked by pronounced tri- 
axial shapes: the average ratio between the smallest and 
largest axis c/a ~ 0.49. T his ratio agrees with w ell with 
the value of 0.45 found by IShandarin et all (|2006h (which 
is perhaps a little surprising given the totally different void 
definitions used in that paper). 

In summary, we find voids to be slightly prolate, with 
average axis ratios in the order of c : b : a ~ 0.5 : 0.7 : I. 

Two important factors contribute to this flattening. 
Even though, inter nally, void s tend to become more spheri- 
cal as they expand (Ickc 1984), perfect sphericity will hardly 
ever be reached. Before voids would be able to become spher- 
ical they would encounter surrounding structures such as 
overdense filaments or planar walls. Even more important 
may be the fact that, for voids, external tidal influences are 
particularly important: voids will always be rather moder- 
ate underdensities since they can never be more underdense 
than 8 — — 1. The external tidal forces drive a significanr 
anisotropy in the development of the voids, and in the ex- 
treme cases may cause complete collapse of the void along 
one or more of its axes. This is another aspect of the way 
in which tidal forces shap e the Cosmic Web, as has been 
emphasized bv lBond et al\ (| 19961 ^) . 



4 VOID ALIGNMENTS 

Large scale influences play a major role not only in shaping 
individual voids but also on their mutual arrangment and 
organization. A particularly important manifestation of this 
is the mutual alignment between voids. In this section wc 
establish the significance of this effect, and in the next sec- 
tion we analyse the correlation between the tidal force field 
and the void orientation. 

4.1 Void Alignments: definitions 

We take the void centers in our sample to define a spa- 
tial point process. Three void alignment measures have 
been investigated. Each i s a marked correlation function 
IIBeisbart fc Kersher|[2OQ0l: IStovan fc Stovanlll994T) and as- 
sesses the degree of alignment as a function of the distance 
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r between the void centres. The mark for each void pair is 
the cosine of the angles between two specific axes of the void 
or a combination of these. 

The first alignment measure is the full shape correlation 
Ass- The second measure concerns the alignment strength 
A33 between the longest axes of the voids and the third 
measure concerns the alignment strength An between the 
shortest axes of the voids. To determine these functions we 
evaluate the alignment between each pair of voids as a func- 
tion of the separatin of their centers r. 

The alignment between a pair of voids m and n, with 
traceless shape tensors S m ,ij and S n ,ij, is specified in terms 
of the following set of quantities, 



T S s(m, n) 

Tn (m,n) 
F22(m, n) 
r 33 (m,n) 



e s i(m) • e a i(n) 
e S 2(m) ■ e a 2(n) 
e s3 (m) ■ e s3 (n) . 



(4) 



in which <S m and S n are the norms of the corresponding 
(traceless) shape-ellipsoids. In essence Fn, 1^2 and are 
the cosines of the angles between the corresponding void 
axes, directed along the shape eigenvectors e a i, e S 2 and e„3. 
These alignment measures do not depend on the degree of 
ellipticity of each void, they focus exclusively on the relative 
orientation of the void axes. 

The corresponding marked correlation functions are an 
expression of the strength of the alignments as a function of 
distance r between the void centers, 

Css(r) = {T S s) (5) 

Cn(r) = (rii), c 2 2(r) = <r 2 2), c 3 3(r) = <r 3 3) 

The brackets indicate the ensemble average over all void 
pairs (m, n) whose mutual separation r lies within a band 
Ar around r, |r m — r n | £ [r,r + Ar]. In practice, we bin the 
void pairs in our void sample into bins which each contain 
5000 pairs. This ensures that each bin has equal statistical 
significance. 

Note that because the angle between the correspond- 
ing void ellipsoid axes of two voids m and n has a value in 
the range 6(m,n) ^ vr/2, the alignment Tu(m,n) be- 
tween these two voids lies within the range ^ Tn(rn,n) = 
cos# mn ^ 1. Each of the axis correlation functions Ca 
(i=l,2,3) is the average cosine between corresponding two 
void axes, Ca = (cos 9). For a perfectly randomly oriented 
sample we would therefore expect Cn = C22 = C33 = 0.5. 
In the case of perfect alignment we would obtain Ca — 1.0. 
A value of Ca « would suggest an anti-alignment, a cross 
pattern, corresponding to an orientation angle 8 mn m tt/2. 

4.2 Void Alignments: significance 

We can gain insight into the significance of the values of the 
alignment measures Ca by comparing the values obtained 
for a sample of vectors in which a perfectly aligned subset 
is mixed with a entirely randomly oriented set of voids. 

In this comparison sample, suppose a fraction f r of the 
void sample has a random orientation while the remaining 
fraction f a = l — f r has a singular orientation along one 



direction. We may then infer the implied random fraction 
f r in the comparison sample for the obtained values for Ca 
in our real void sample. 

It is straightforward to infer the implied values of the 
alignment measures Ca in the comparison sample: 



Ca = (cosfl) 



= 0.5/ r 2 + frfa + fl 

= 0.5 if r - l) 2 + 0.5. 



(6) 



For example, if we were to have 20% aligned voids mixed 
with 80% random ones, the alignment would be Ca = 0.52. 
An alignment of Ca = 0.63 would correspond to « 50% 
aligned and « 50% random voids, while Ca = 0.75 would 
imply « 70% aligned and « 30% random ones. 

4.3 Void Alignments: results 

Figure presents the results for the measured correlations 
of the void orientations. Each curve corresponds to the ori- 
entation correlation between one of the three void axes of 
voids: the solid black curve to that of the longest axis, C33, 
the dot-dashed curve to that of the shortest axis, Cn, and 
the intermediate axis C22 (dashed curve). 

The figure shows the longest and the smallest axis ex- 
hibit the strongest alignment. At short distances their align- 
ment is very strong and remains substantial out to rela- 
tively large radii of 20 — 30h~ Mpc. Beyond this distance 
the alignment rapidly declines towards random values. At 
the smallest distances probed the alignment of the longest 
and shortest axis reaches values of around 0.65. This corre- 
sponds to a random component of only 30 percent, in combi- 
nation with 70 percent perfectly aligned voids. By contrast, 
the second axis shows a striking lack of alignment on all 
scales and, except at the shortest separations, appears to be 
mostly randomly oriented. This is partially due to the con- 
siderable scatter in the orientation of the intermediate axis. 
This may be a reflection of the crudeness of approximating 
the void shapes by ellipsoids. However, even taking account 
of the implied scatter the intermediate axis remains more 
weakly correlated. 

To provide an estimate of the significance of the results 
we have estimated the standard deviation expected for a 
sample of perfect randomly oriented voids. The 2a spread 
is depicted by means of the coloured bar around the expec- 
tation value for a random sample, Ca = 0.5. We obtained 
these estimates by a Monte Carlo experiment in which we 
randomly shuffled the shape ellipsoids over the the voids 
and subsequently re-measured the alignment. Note that the 
variance remains nearly constant over the whole range. This 
is a consequence of our decision to define bins having equal 
number of particles per bin. 



5 THE CAUSE OF VOID ALIGNMENT 

The possibility that the alignment may be an artefact of the 
WVF method has to be considered: large spherical voids 
that were chopped in half would provide an anisotropic void 
set in which neighbouring voids were strongly aligned. How- 
ever, this it can quickly be eliminated simply by examining 
figure \T\ which shows no such phenomenon. We shall later 
on present further evidence based on tidal fields that WVF 
is not itself the source of the alignment. 
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Figure 3. The alignment of the three void axes. The figure shows 
the (marked) correlation functions Cn(r) (solid), C22(r) (dashed) 
and C33(r) (dot-dashed), the ensemble average of the cosine of the 
angle between the longest, medium and shortest axes of voids at a 
distance r. A pure random orientation corresponds to Ca(r) = 0, 
a perfectly aligned one Cu{r) = 1. The coloured bar around the 
value Cu = 0.50 indicates the 2 — a dispersion of a similar random 
sample. 



The reason for the alignment between voids may be due 
to any of a number of physical effects. Here, we consider 
three possible physical sources for this alignment: 1) the 
initial conditions, 2) the geometric packing of the voids and 
3) the large scale tidal field. 

5.1 Initial conditions 

The first possibility, that the structural correlations are born 
with the primordial density field, is difficult to address since, 
at least in the linear regime, today's matter distribution a 
direct reflection of those initial conditions. Thus the density 
fluctuations and the tidal fields that distort them are both 
present initially, either being deducible from the other. The 
Cosmic Web, voids, galaxy halos as well as the tidal field are 
all manifestations of the same initial matter distribution. 

Voids emerge out of a primordial Gaussian density field. 
Density peaks a nd dips in a Gaussia n random field are not 
only cluster ed ( Bardeen et al.l Il986t ) , they are also mutu- 
ally aligned (|Bondlll987h ~ This alignment stretches out over 
scales over which they possess a nonzero correlation func- 
tion, approximately the same range as that of the compa- 
rable population of rich clusters. Even when their evolution 
would not be influenced by external forces the voids would 
retain the memory of their initial configuration. Their ini- 
tial shape would be (slightly) attenuated by their internally 
driven expansion and the void would keep its initial orien- 
tation. 



The primordial alignment of density peaks is enhanced 
by subsequent nonlinear interactions with surrounding mat- 
ter condensations. Voids, however, will be far less affected 
due their relatively large spatial extent and relaively low 
density contrast. If anything, their initial alignment will be 
strengthened by the influence of the large scale tidal field. 

5.2 Void packing 

Voids are a unique constituent of the cosmic inven- 
tory in that they represent a volume-filling population 
jSheth fc van de Weygaert! 12004 ). Consequently, the place- 
ment of voids in the large scale matter distribution is 
a nontrivial issue related to the random object packing 
problem in mathemati cs and material's science (see e.g. 
IConwav fc Sloandl 19981 ). 

We have argued before that an important consequence 
of the tight packing of voids concerns their shape. As they 
try to expand outward they will meet up with the surround- 
ing large scale structures. This involves the surrounding fila- 
ments, walls and clusters as well as their neigbouring voids. 
This strongly restrict s their tendency to become mor e spher- 
ical as they expand (|lckd|l984l; iBertschingerl ll985T). While 
the i nner cores of voids (|van de Weygaert fc van KamperJ 
1993) remain spherical, the surroundings mould them into 
a generally more flattened configurations. (The shape of re- 
gions defined by various level-sets in the i nitial density field 
is discussed in detail bv lShandarin et all (|2006l )). 

Mathematical studies of the packing problem focus 
mainly on the packing of equivalent spheres. The more com- 
plex problem of ellipsoid packings was recently adressed 
in a seminal experimenta l study on the tig ht packing of 
(oblately shaped) M&M's ijDonev et al.ll2004i ). Not only did 
they prove that such ellipsoids allow a tighter packing than 
feasible with purely spherical objects, but also that this goes 
along with a configuration of strongly aligned ellipsoids (see 
fig. lB. lDonev et~aill2004 ). 

The situation for the packing of voids is considerably 
more complex. Cosmic voids have, as we have shown, a 
considerable diversity in size and shape. Moreover, unlike 
M&M's, voids are not solid bodies: the void ellipsoids may 
not be entirely contained within the boundaries of their void 
patches and as a result may overlap. Nonetheless, we might 
on purely geometric grounds expect an alignment of a void 
with its "natural" neighbours and that, locally, the orienta- 
tion of voids is strongly constrained by the shape and orien- 
tation of its neighbours. 

In order to further assess the situation we have analyzed 
four different configurations: a ACDM initial Gaussian field, 
a Poisson field, a regular crystalline configuration and an 
irregular Voronoi cellular configuration. These reflect one or 
more aspects of the genuine void distribution (an accompa- 
nying publication will provide the details of the study). 

In each of these special configurations we see a rapid 
decline in alignment strength beyond the characteristic void 
size Ry, even turning into anticorrelation around 2Rv- 

Figure[3]shows the large scale correlations for void align- 
ment in the N-Body model we have studied. We have marked 
the scales Rv and 2Ry by means of two vertical lines. The 
fact that in the void sample of this study we find signifi- 
cant alignments beyond 2Rv is a strong argument for the 
alignments being due to other factors. 
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5.3 Large Scale Tidal Forces 

The intricate patterns that constitute the cos mic web are 
a dir ect product of the large scale tidal field (jBond et al\ 
1996). Voids are an integral part of this pattern and their 
shape and orientation is expected to react strongly to the 
surrounding inhomogeneous matter distribution via their 
tidal influence. The tidal field is a reflection of the underlying 
inhomogeneous cosmic matter distribution and contains con- 
tributions from matter inhomogeneities over a large range of 
scales. 

In a sense, the influence of the tidal field on void ori- 
entations is strongly coupled to the initial conditions (see 
above): already in the initial density field the local shape 
of a density trough is considerably correlated with the tidal 
field configurations over a large range of scales. 

The key issue is the influence of tidal forces on the evo- 
lution of the shape and orientation of voids. This includes 
the impact on the nonlinear development of the void. We are 
particularly interested in identifying the largest scales over 
which we can recognize a clear alignment between the tidal 
field and the orientation of the voids. Their spatial coherence 
will be crucial for the void- void alignments. 

In the next section we will focus on this specific issue. 



6 INFLUENCE OF TIDAL FIELD 

To test the relationship between the large scale tidal force 
field and the void correlations of which they are supposed to 
be the source we investigate the correlation, at the location 
of each of the voids, between the tidal tensor %j and the 
shape tensor iSy. 



6.1 Tidal Field 

We compute the traceless tidal shear components, 



+ 1 1 



dridri 



(7) 



for the entire simulation volume. The gravitational potential 
4>, due to the underlying density distribution, is determined 
from the DTFE density field 5(r) by solving the Poisson 
equation in Fourier space. 

The eigenvalues of the tidal tensor are ordered as ti > 
t2 > £3, with corresponding unit eigenvectors e t i, e t 2 and 
e t 3. We also draw a distinction bewteen the compressional 
component of the tidal field, that part of the tidal field corre- 
sponding to the positive eigenvalues of 7y , and the dilational 
component. The compressional component always includes 
the component tieti while the dilational component always 
includes i3et3. Depending on its sign, t 2 et2 may contribute 
to either the compression or to the stretching of the mass 
clement. 



6.1.1 Tidal-Tidal alignments 

In our bid to identify the source of the void-void shape cor- 
relations we also need to understand the auto-correlations 
between the tidal field contributions originating at different 
scales. To this end we define a tide-tide alignment measure 
Att{Ri,R2). 



At a location r in the density field we may compare the 
orientation of the tidal field filtered at two different scales 
Ri and R2, yielding the alignment 



T T T(r;Ri,R 2 ) = 



]T T ij {r;R 1 )T ij {T,R 2 ) 

— % i3 

T(r; R±) T{r;R 2 ) 



(8) 



where T(v;R\) is the norm of the tidal tensor %j{v;R\) 
filtered on a scale R\ and T(r; R2) on scale i?2- 

The ensemble average of Ttt(t; Ri, R2) is the tide-tide 
alignment ^^(^1,^2), 



Att(Ri,R 2 ) = (T T t(Ri,R 2 )) 



(9) 



This quantity is computed by taking the average of 
TTT(r; Ri, R 2 ) over a large number of locations r. 



6.1.2 Void-Tidal Field alignment: formalism, 

In order to trace the contributions of the various scales to 
the void correlations we investigate the alignment between 
the void shape and the tidal field smoothed over a range of 
scales R. The smoothing is done in Fourier space using a 
Gaussian window function W*(k;R): 



%j(r;R) 
-ilH 2 



(10) 



dk 

(2^ 



fc-' 



W*(k;R) <5(k) e - ik 



Here, <5(k) is the Fourier amplitude of the relative density 
fluctuation field at wavenember k. 

We evaluate the alignment Ats(R) between the void 
shape ellipsoid and the local tidal field tensor Tij(R), Gaus- 
sian filtered on a scale R at the void centers. This is done by 
evaluationg, for every void, the function IYs (m, R) at the 
void centers: 



r TS {m;R) = 



T ZJ {r m ;R) 



(11) 



S m T(r m ;R) 

where T(r m ; 7?) is the norm of the tidal tensor 7ij(r m ) fil- 
tered on a scale R. The void-tidal alignment Ats(R) at a 
scale R is then the ensemble average 



Ats(R) 



(^ts(R))- 



(12) 



which we determine simply by averaging rYs(m, R) over the 
complete sample of voids. 



6.2 Void-Tidal Field alignment: results 

A visual impression of the strong relation between the void's 
shape and orientation and the tidal field is presented in the 
lower righthand panel of fig. [1] The tidal field configuration 
is depicted by means of (red-coloured) tidal bars. These bars 
represent the compressional component of the tidal force 
field in the slice plane, and have a size proportional to its 
strength and are directed along the corresponding tidal axis. 
The bars are superimposed on the pattern of black solid 
watershed void boundaries, whose orientation is emphasized 
by means of a bar directed along the projection of their main 
axis. 

The compressional tidal forces tend to be directed per- 
pendicular to the main axis of the void. This is most clearly 
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Figure 4. Solid line: the local alignment Ats(R) of voids with 
the tidal field, as a function of the smoothing radius R of the tidal 
field. Dashed line: alignment Alp^(R s ,R) between tidal field on 
scale R s = 6h~ 1 Mpc and tidal field on scale R (see text). The 
relation is based on the prediction for a Gaussian random field 
with a ACDM spectrum. Dot-dashed line: tide-tide alignment 
Att {Us, R) measured from the N-body simulation described in 
the text. 



in regions where the forces are strongest and most coher- 
ent. In the vicinity of great clusters the voids point towards 
these mass concentrations, stretched by the cluster tides. 
The voids that line up along filamentary structures, marked 
by coherent tidal forces along their ridge, are mostly ori- 
ented along the filament axis and perpendicular to the local 
tidal compression in these region. The alignment of small 
voids along the diagonal running from the upper left to the 
bottom right is particularly striking. 

To quantify and trace the tidal origin of the alignment 
we have plotted in figure U the local shape-tide alignment 
function Fts versus the smoothing radius R. By smoothing 
over a scale R we suppress the contribution by mass concen- 
trations and structures with a scale smaller than R. Like- 
wise, we have evaluated the tide-tide alignment Att(Ri, R) 
as a function of the same variable filtering scale R. The fil- 
ter scale R\ (eqn. [9} is fixed at a value of Ri ~ 6ft" 1 Mpc, 
slightly smaller than the mean void size. At this scale we 
expect a strong influence of the local shape on the tidal 
field, so that this quantity is expected to provide additional 
information on the dynamical influence of the large-scale 
tidal field on the dynamics of the void. Note that this im- 
plies there to be a slight offset between the actual mean 
void size, R v « 8.5/i -1 Mpc and the peak scale of Ats- The 
latter is biased towards the more strongly clustered small 
voids, whose alignment is often due to one coherent nearby 
large-scale structure, often a filament or cluster. 



The following observations can be made: 

• The alignment is strong over the whole range of smoothing 
radii out to R « 20 - SO/t^Mpc (Ats ~ 0.25). 

• The alignment peaks at a scale of R ~ 6/i _1 Mpc. This 
scale is very close to the average void size (see fig. f2J), 
and also close to the scale of nonlinearity. The latter is 
not a coincidence: the identifiable voids probe the linear- 
nonlinear transition scale. 

• There is still a remarkably strong alignment signal at radii 
larger than R > 20/i _1 Mpc (where Ats ^ 0.3). It forms a 
strong argument for the substantial role of the large scale 
tidal field in aligning the voids. If there were no external 
large-scale contribution the alignment Tts would have 
disappeared for smoothing radii R in excess of the average 
void scale. 

• The fact that the alignment Fts at the peak scale R ~ 
6/i _1 Mpc is less than twice that of at R ~ 25h~ Mpc 
shows that local contributions, in particular due to the 
voids' own ellipticity, are only partially responsible for 
the shape-tidal field alignment. 

We have also analyzed higher resolution simulations. The 
results we found here carry over to the higher resolved sim- 
ulations. The smaller box of the latter, however, contain 
fewer voids and there is a smaller range over which the in- 
fluence of tidal fields can be included. The results do remain 
significant but are marked by more noise. 



6.3 Tidal Connections 

We have established a strong correlation between void ori- 
entation and the large scale tidal field. We now compare 
this with the self-alignment of the large scale tidal field that 
would be expected on the basis of linear perturbation theory. 
The dashed curve in fig.|4]is the linear theory tide-tide align- 
ment A T t{Rs,R), with R s = 6/i~ 1 Mpc. On scales greater 
than the void scale the void-tidal field alignment follows the 
tidal field alignment. On smaller scales, however, we start 
to see a marked deviation: the alignments are stronger than 
expected from the initial linear field. This may be because 
of nonlinear effects on small scales: the tidal field on the 
smaller scales where the density fluctuaions are nonlinear is 
even strong enough to cause the destruction of small voids 
near the filaments of the c osmic web where there is a st rong 
anisotropic force field (see Ivan de Weygaert et aZ.|[20oi ). 

We also measured the tide-tide alignment Att (Rv , R) 
for the tidal field of our ACDM cosmology where any non- 
linear effects that are present are accounted for. The dot- 
dashed curve is the result (note that the amplitude has been 
scaled so that the amplitudes of both Att curves match the 
peak of the Ats curve). The fact that, on small scales, there 
is clear agreement between the void-tide alignment and the 
tide-tide alignment underlines the important role of tidal 
forces throughout the formation history of structures in the 
cosmic web. 

As an incidental remark, the existence of the long range 
tide - void alignment tells us that the void anisotropy is not a 
mere artefact of the WVF method chopping spherical voids 
in two (which would of course lead to non-spherical voids 
having correlated alignments). 
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7 CONCLUSION AND DISCUSSION 

We have investigated the shapes and alignments of voids 
found in the large scale matter distribution of the VLS GIF 
N-body simulation of structure formation in a ACDM uni- 
verse. The void sample has been identified using the Water- 
shed Void Finder (WVF) technique. 

7.1 Overview of Results 

Although voids tend to be less flattened or elongated than 
the halos in the dark matter distribution, they are never- 
theless quite nonspherical: they are slightly prolate with 
axis ratios on the order of c : b : a ~ 0.5 : 0.7 : 1. Two 
important factors contribute to this flattening. Inter nally, 
voids tend to become more spherical as they expand l|lckel 
Il984f ). However, they will never be able to reach perfect 
sphericity before meeting up with surrounding structures. 
Furthermore, external tidal forces will contribute substan- 
tially to the anisotropic development of the voids. In this 
study we have established the influence of the tidal fields 
as being an important agent in the dynamical evolution of 
voids. Indeed, our study provides strong confirmation that 
tidal forces dominate t he shaping of the Cosmic Web, as 
emphasized long ago bv lBond et all (|l996T l. 

We have also investigated the relative orientation of the 
voids. The orientation of voids appears to be strongly cor- 
related with alignments spanning distances > 30/i _1 Mpc. 
This is true of both the shortest and longest axes of the 
void shape. This coherence is quite conspicuous in plots of 
the density distribution. To test for such void alignments in 
observational data a follow-up study will apply the WVF 
formalism to galaxy redshift surveys. 

We find an intimate link between the cosmic tidal field 
and the void orientations. Over a range of scales we find 
a strong alignment of the voids with the tidal field aris- 
ing from the smoothed density distribution. Given that the 
alignment correlations remain significant on scales consider- 
ably exceeding twice the typical void size, our results show 
that the long range external field is responsible for the align- 
ment of the voids. This confirms the view that the large scale 
tidal force field is the main agent for the spatial organization 
of the Cosmic Web. 

7.2 Insights 

We have argued that the large scale tidal field is not only 
the main agent responsible for the shapes of the voids but 
also for the resulting alignments of the voids. Locally, the 
orientation of a void turns out to be strongly aligned with 
the tidal force field generated by structures on scales up to 
at least 20 — 30/i~ 1 Mpc. On scales comparable to the average 
void size we also found indications of nonlinear effects, with 
the void's orientation reacting to the small-scale nonlinear 
influences. 

This conclusion agrees wi t h tha t of a similar study of 
halo alignments by iLee et all (|2007h . However, the align- 
ment of halos tends to be strongly attenuated by local non- 
linear effects (e.g. van Haarlem fc van de Weygaertl 1 19931 ; 
lAragon-Calvo et al.l 20071 ) rendering the final halo alignment 
signal weaker than that for voids. 

One aspect we have not yet adressed in detail is the 



entanglement of initial conditions and the influence of the 
tidal field. On large scales, where the density fluctuations 
are in the linear regime, these are different aspects of the in- 
homogeneous matter distribution: if one is known the other 
can be determined. In this respect the indication of nonlin- 
ear effects in the tidal- void alignment Ats (fig- H]) at small 
scales is interesting. They may very well be an expression 
of the strong influence of coherent overdense filaments and 
dense compact clusters on the evolution of small voids in the 
out er regions of large voids . 

Ivan de Weygaert et all l|2004l 'l argued that the col- 
lapse of small voids, an essential aspect of the hierar - 
chical evolution of voids (|Sheth fc van de Weygaertl I2004T I. 
manifests itself in a tidally induced anisotropic contrac- 
tion of small underdensities at the void boundaries. The 
dynamics and the fate of these small underdensities is 
often decisively influenced by external anisotropic forces 
(van de Weygaert fc Babul|[l996l ). While such forces may af- 
fect the nonlin ear evolution of halos and clusters to a c onsid- 
erable extent l|Bond fc Mversll 19961 ; Isheth et al.ll200lf ). they 
do not go as far as deciding their fate. 

This work is part of a study of the evolution of voids un- 
der the influence of external influences. A proper framework 
for following the gradual hierarchical buildup of void regions 
under the influence of ext ernal tidal forces is provided by the 
void-patch description (|Platen. van de Weygaert fc Jonea 
120081 ; iBond et aZlll996T ). To some extent the void-patch for- 
malism is a more accurate approximation of reality than 
for the equivalent peak-patch: because of their approximate 
uniformity void evolution is accurately described by the el- 
lipsoidal model. 

A perhaps even more complex issue is that of the iden- 
tity of the outer regions of voids, the regions where the void 
merges into the surrounding matter distribution. We have 
argued that this will play an important role in determin- 
ing the shape of voids. However, a proper definition of a 
void boundary does not (yet) exist. Arguably the WVF has 
proven to represent a major advance along these lines and 
we should soon be able to be more definite on the issue. 

Having established the external influences on a void's 
evolution, we may try to understand the sensitivity of the 
void population to the underlying cosmology. This depends 
to a considerable degree on th e differences in sc ale depen- 
dence of the tidal force fields. ILee fc Parkl J2007) did find a 
considerable influence. This will open the door to the use of 
voids in determining the global cosmology. 
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